LINEAR STABILITY ANALYSIS FOR TRAVELING WAVES OF SECOND 

ORDER IN TIME PDE'S 
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Abstract. We study traveling waves ip c of second order in time PDE's uu + Cu + N(u) = 0. 
The linear stability analysis for these models is reduced to the question for stability of quadratic 
pencils in the form X 2 Id + 2c\d x + H c , where H c = c 2 d xx + C + N'(ip c ). 

If H c is a self-adjoint operator, with a simple negative eigenvalue and a simple eigenvalue at 
zero, then we completely characterize the linear stability of cp c . More precisely, we introduce 
an explicitly computable index lo*(H c ) G (0, oo], so that the wave ip c is stable if and only if 
|c| > ui*(T-L c ). The results are applicable both in the periodic case and in the whole line case. 

The method of proof involves a delicate analysis of a function Q, associated with W, whose 
positive zeros are exactly the positive (unstable) eigenvalues of the pencil \ 2 Id + 2c\d x + H. We 
would like to emphasize that the function Q is not the Evans function for the problem, but rather 
a new object that we define herein, which fits the situation rather well. 

As an application, we consider three classical models - the "good" Boussinesq equation, the 
Klein-Gordon-Zakharov (KGZ) system and the fourth order beam equation. In the whole line 
case, for the Boussinesq case and the KGZ system (and as a direct application of the main 
results), we compute explicitly the set of speeds which give rise to linearly stable traveling waves 
(and for all powers of p in the case of Boussinesq) . This result is new for the KGZ system, while 
it generalizes the results of [2] and [1], which apply to the case p — 2. For the beam equation, 
we provide an explicit formula (depending of the function H^cIL 2 ); which works for all p and for 
both the periodic and the whole line cases. 

Our results complement (and exactly match, whenever they exist) the results of a long line of 
investigation regarding the related notion of orbital stability of the same waves. Informally, we 
have found that in all the examples that we have looked at, our theory goes through, whenever the 
Grillakis-Shatah-Strauss theory applies. We believe that the results in this paper (or a variation 
thereof) will enable the linear stability analysis as well as asymptotic stability analysis for most 
models in the form u u + Cu + N(u) — 0. 



1. Introduction and motivation 
The main motivation of our study is the following abstract second order in time nonlinear PDE 

(1) utt + C x u + N(u) = 0, (t,x) GR+ xR d or (t, x) G R 1 x [-L, L] d , 

where C x is a given linear operator, acting on the x variable and N(u) is the nonlinear term. We 
outline some relevant examples of interest, see Section 1.1 below, but notice that we do consider 
both periodic boundary conditions as well as vanishing at infinity solutions. These two scenarios 
will be considered simultaneously, since our method works equally well in both cases. 
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Our interest is in the study of the stability properties of traveling waves, that is, solutions in 
the form ip{x + ct). Clearly, these satisfy the stationary PDE 

d 

(2) C x ip + ^ CiCj'PxiVxj + N((p) = 

We take the ansatz u = tp(x + ct) + v(t,x + ct) and plug it in (1). Taking into account (2) and 
dropping all quadratic and higher order terms in v, we arrive at the following 

d 

(3) vtt + 2 (c, V x v t ) + c jCk | + Cv + N'(ip)v = 

j,k=l ax 3° X k 

Therefore, introducing the operator Hg = C x + Ylj k=i c j c kd Xj d Xk + N'(ip), we are lead to study 
the following problem 

(4) v tt + 2 (c, V x v t ) + Hu = 0. 

This is the linearized stability problem for the nonlinear equation (1) in a vicinity of the special 
solution (p(x + ct). Before we move on with our exposition, we provide some informal definitions 
and discussions to motivate our interest in the problem. There are several notions of stability that 
are of interest. Generally, linear stability is easier to check than (and is often a prerequisite to) 
nonlinear stability. Even within the linear stability, we distinguish between spectral and linear 
stability. Namely, for an evolutionary problem in the form zt = Mz, where M is a closed operator 
generating Cq semigroup, we say that we have spectral stability, if a(M) C Z- = {A : KA < 0}. 
We say that the same problem has linear stability, if the solutions grow at infinity slower than any 
exponential 1 , i.e. for every 5 > 0, lim^oo e _5 *||z(t)|| = 0. The relationship between spectral and 
linear stability has been well-explored and documented in the literature, and we will not dwell 
on it, except to point out that in principle (and in the presence of the so-called spectral mapping 
theorem for the generator M), these are equivalent and amount to lack of exponentially growing 
modes, that is solutions in the form Aitp = Xip for A : 3ftA > 0. 

On the other hand, we have two distinct notions of nonlinear stability - orbital and asymptotic 
stability. Assuming for simplicity that the only invariance of the system is translation, orbital 
stability requires that a solution for the full nonlinear equation that starts close to a traveling 
wave stay close for all times to a (time-dependent) translate of the starting wave. Asymptotic 
stability requires a bit more, namely that the perturbed profile will actually converge to a (time- 
dependent) translate of the wave. 

There is a large body of literature that deals with this problem in various models. We would 
like to point out that the powerful methods of Grillakis-Shatah-Strauss, [14], [15] reduce (in most 
cases) the problem of orbital stability to checking certain conditions on the linearized functionals. 
We mention in this regard the papers [4], [7], [25], [26], which treat models of interest to us in 
this paper. We also note that establishing orbital instability for a given problem seems to be 
harder and requires more problem specific efforts, [25]. 

Our interest here is the question of linear stability of such traveling waves <p c , that is, whether 
there are exponentially growing solutions of (4) in the form e xt ip(x). We pose the following 

Question: For which values of c, the corresponding traveling wave ip c determined by (2) is 
linearly /spectrally stable? More precisely, for which c, the equation (4) has a solution in the form 



Here, note that one cannot require bounded orbits. Indeed, even in the ODE case, if one takes a Jordan block 
A of size d, e tA has polynomial growth t d 
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Note that this question has been studied thoroughly in the last twenty years. Here, the methods 
are completely different, although the results in the end must be related with those obtained by 
orbital stability methods. After all one expects orbitaly stable waves to be linearly stable and 
vice versa (although this sometimes fails at points where the stability character changes). 

The Evans function method has been used in connection to the linear stability of such waves, 
[2], [30], [31], [23]. However, the method has well-known limitations, in particular dealing with 
systems, which is the issue at hand here, since the equations are second order in time. In fact, 
the papers quoted here are among the few that deal with second order in time equations. 

Another method that has been developed is the method of "indices counting" , which basically 
relates the number of unstable eigenvalues of self-adjoint entries of M. to the number of unstable 
eigenvalues of M itself. This has been mostly useful in the KdV and Schrodinger type systems, 
but it has also played important role in general Hamiltonian and dissipative systems, [17], [18], 
[19], [9], [10]. We would like to point out that some of these results have enabled the consideration 
of spatially periodic waves, [10], which is one of the goals of this project as well. 

A third method for establishing mostly sufficient conditions for instability is through a direct 
construction of unstable modes, [25], [24]. In fact, there are many other impressive results 
in the literature, mostly for standing waves, which provide "instability by blow up" for such 
waves. These are done typically by constructing clever Lyapunov functionals, associated with the 
unstable modes. 

Our goal in this paper is to develop a fairly general and systematic theory, which treats the 
question for linear stability (i.e. the existence/nonexistence of such -ip) of traveling waves of 
second order in time systems. Indeed, we present a complete answer to this question, in the case 
of (4), when H is self-adjoint, with at most one negative eigenvalue and d = 1. 

1.1. Examples. We consider three basic examples that fit this category, although there are 
numerous others to which our results will be applicable. The reader should also bear in mind 
that one might consider both periodic and whole line boundary conditions, for the models that 
we list below. To keep the discussion simple, we will mostly stick to the whole line case, the 
periodic cases will be addressed in a subsequent publication, [16]. 

1.1.1. The "good" Boussinesq models. Our first example is the "good' Boussinesq-type model, 

(5) u tt + u xxxx - u xx + {u p ) xx = 0, (t, x) G R+ x R 1 

This model was considered by Bona-Sachs, [4] as a model for propagation of small amplitude, 
long waves on the surface of water. This is indeed an equation that belongs to a family of 
Boussinesq models, which all have the same level of formal validity, however (5) exhibits some 
desirable features, like local well-posedness in various function spaces, [4]. Interestingly, global 
well-posedness for (10) does not hold, even if one requires smooth initial data with compact 
support. In fact, there are "instability by blow-up" results for such unstable traveling waves for 
this equation. 

It is easy to see that there exists one-parameter family of traveling waves of the form ip(x — 
ct), \ c\ £ (—1, 1), which obey the equation 

(6) c 2 ip + ip" - if + <p p = 
and which have the explicit form 



P+1 ^(l-c< 



P-1 2 / s/l - c 2 (p- 
sechp- 1 I 
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One should of course recognize the standard sech solitons of the generalized KdV hierarchy. 
This is not surprising, since the governing ODE for the traveling waves of (10) is related to the 
corresponding ODE for gKdV. 

In [4], the authors proved orbital stability for these solutions, provided p < 5, ^ p 2 1 < \c\ < 1, 
while nonlinear instability was established in [25]. Both of these results were achieved via 
the methods of the Grillakis-Shatah-Strauss theory (developed in [14], [15]) for orbital stabil- 
ity/instability. 

1.1.2. The Klein- Gordon- Zakharov system. Consider the following system of coupled wave equa- 
tions 2 



(7) 



utt — Uxx + u + nu = (i, x) G R+ x R 1 
n tt ~ n xx - \{\u\ 2 ) xx = 0, 

which describes the interaction of a Langmuir wave and an ion acoustic wave in a plasma 3 , [13]. 
Taking the traveling wave ansatz u(t,x) = ip(x — ct),n(t,x) = ip(x — ct), we derive from the 

second equation that (c 2 — l)?/>" = ^(p 2 )", whence (c 2 — = since the functions p>, if) decay 
at infinity. Using this relation in the first equation yields the following ODE for p 

(8) - (1 - cV ' + ^-2Tl^) =0 

It is well-known (and in fact can be checked easily, based on simple rescaling arguments in the 
ODE's governing traveling wave solutions for the KdV equation) that (7) admits an one parameter 
family of traveling wave solutions in the form u(t, x) = (p(x—ct),n(t, x) = tp(x-ct) for c G (—1, 1), 
where 



(9) 



<p(y) = 2\/l — c 2 sech 

= -2sech 2 ( 7f L=) 



In fact, there exists a two-parameter system of solitary-traveling waves of which the family dis- 
played above is a particular case, [7]. In the same paper, the author shows orbital stability for 
such waves, provided 1 > |c| > Below, we solve completely the question for linear stabil- 
ity. Namely, for \c\ G [-^,1) we have linear stability, whereas for |c| G [0, ^) we have linear 
instability. 

1.1.3. The nonlinear beam equation. Another relevant example to consider is the so-called beam 
equation 

(10) u tt + A 2 u + u-\u\ p - 1 u = 0, (t,x) G R 1 x H d or (t,x) G R 1 x [—L,L] d , 

where p > 1, L > and we either require periodic boundary conditions (in the case x G [— L, L]) 
or vanishing at infinity for x G R rf . 

This equation has been studied extensively in the literature. It seems that the earliest work 
on the subject goes back to [32], where (10) is proposed as a model of a suspension bridge. In 
this paper, the authors have also showed the existence of some traveling wave solutions, that is, 
solutions in the form ip(x + ct). 



2 Here the coefficient | in front of the (ImI 2 )^^ can be taken to be arbitrary by rescaling, but the particular choice 
of | will be convenient in the spectral analysis 

3 namely, the complex u describes the fast scale component of the electric field, while the real valued n stands 
for the deviation of ion density. 
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In [26] it was proved that traveling wave solutions exist, in the whole space context, whenever 4 
|c| G (0, y/2). The proof is variational in nature, and this allowed the author to derive conditions 
for orbital stability /instability of such waves for different speeds. In short, the conclusion was that 
such traveling waves are orbitally unstable for small speeds, while an orbitally stable solutions 
are observed for values of the parameter | c| ~ \^2, | cj < y/2. Various other works have explored 
the decay and scattering properties of the linear beam equation uu + A 2 + u = 0, [27, 28, 29]. 

1.2. Some general features of the operators H. In this section, we construct the operators 
H that arise in (4) and discuss some of their general features. Of particular importance will be 
the fact that they are self-adjoint and the fact that they have exactly one negative eigenvalue. 

1.2.1. The beam equation - spectral picture. Clearly, such a function <p will satisfy the ODE 

d 

(11) AV+ c i c Wx iXj +¥>- M P "V = 0, x£(-L,L) d , 

which is supplied by the usual periodic boundary conditions at —L, L (or vanishing at infinity for 
the whole line case). Note that we do not necessarily look for positive solutions of (11), nor do 
we expect to find any in general (because of the lack of maximum principle for A 2 , the ground 
state are not positive) , thus the presence of the absolute value in the non-linear term. In the one 
dimensional case, which will be the main subject of our investigation, (11) reduces to 

(12) c 2 p" + p"" + p-\p\ p ~ l p = 0, -L<x<L. 

Clearly, the question for the existence of such 2L periodic solutions ip is a non-trivial one, but we 
refer the reader to the upcoming publication [11], where it will be addressed. Note though that if 
p is constructed in the space L 2 er [— L,L], then one can automatically bootstrap its smoothness 
to infinity, thanks to the form of (12). 

Even less obvious issue is the linear stability of such waves. In order to simplify matters, 
let us assume henceforth that d = 1 and p is an odd integer (in particular, one should keep in 
mind the model case p = 3). Now, for any such periodic wave <p satisfying (12), set u(t,x) = 
ip(x + ct) + v(t, x + ct), which we plug in (10). After ignoring all the quadratic terms in v, taking 
into account (12) and changing variables x — ct — > x, we obtain the following linearized equation 
around the traveling wave profile (p 

(13) v tt + 2cv tx + v xxxx + c 2 v xx + v - pp p ~ l v = 0, -L < x < L. 
This is then a proper time to introduce the self-adjoint operator H 

(14) H:=d£ + <?% + l-p<f-\ 

with the domain D(%) := H* er [— L, L\. Note that T-L = U(c,p, p) depends implicitly on the speed 
c, the parameter p and the wave p> itself. 

Another fact that should be stated at this point is that % has zero on its spectrum, with 
eigenfunction p' . Indeed, taking a spatial derivative in (12) (under the assumption that p is an 
odd integer) implies H[p'] = 0. Finally, the bottom of the spectrum of T-L is a negative number 5 . 
Indeed, 

(Hp,p) = -(p-l)J p p+l (x)dx<0, 



as we shall see in our arguments later on, this is a natural space of parameters 
'note that Weyl's theorem, in the whole line case <J a . c .(H) = [l,oo) 
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hence <j('H)n(— oo, 0) / 0. Unfortunately, we cannot verify that this eigenvalue is simple, although 
this is easily seen in numerical simulations, for all values of the parameter c. 

Note that the computations above for the beam equation in the periodic case apply equally 
well for the whole line homoclinic solution produced by Levandosky, [26]. 

1.2.2. The "good" Boussinesq model: spectral picture. Here we outline the relevant operators for 
the Boussinesq model. As we shall see though, the scheme outlined in (4) does not quite work 
and one needs a little trick to make it work. 

In any case, if we impose the ansatz u = (p c (x + ct) + v(t,x + ct) and we get the equation 
vu + 2cvtx + Tv = 0, where 

Tv = & x v - (1 - c 2 )d 2 x v +p(<pe- 1 v) xx 

Note that the operator T is not self-adjoint and thus unsuitable for our method. However, if we 
introduce the variable z : z x = v, we get the following linearized equation in terms of z, 

(15) z ttx + 2cz txx + T[z x ] = 0. 
Note that T[z x ] = d x [H[z]], where 

(16) Hz = d x z - (1 - c 2 )d 2 x z + p^- x z x ) x 

Thus, the linearized equation becomes d x [zu + 2czt x + Hz] = 0, which in view of the conditions 
lim^i^^ |<9°z(x)| = 0, implies z tt + 2cz tx + Hz = 0. Thus, the linearized equation to consider is 
again z tt + 2cz tx + Hz = 0, where H is defined in (16). As we prove later on, H has one simple 
eigenvalue at zero and one simple negative eigenvalue. 

Note that while it is clear the presence of unstable mode for zu + 2czt x + Hz = implies linear 
instability for (15) (via v = z x ), the reverse it is not immediately clear. The reason for that is 
that the operator d x is not invertible and one may potentially have a solution e xt V{x) of (15), 
while the corresponding "solution" e xt Z{x) = e xt d~ 1 V{x) may not be an L 2 function. 

2. Main Results 

In this section, we present our results. We start first with a abstract form of the linearized 
problem, which involves the theory of quadratic pencils. 

2.1. Stability /instability results for quadratic pencils. In connection with (13), we may 
write it schematically in the form 

v tt + 2cv tx + Hv = 0. 

Definition 1. We say that the periodic wave (p is linearly unstable, if there exists A : > ; 
and a function ip, so that then the following equation is satisfied 

(17) \ 2 ip + 2c\ip x + Uip = 
Otherwise, we say that the traveling wave <p is linearly stable. 

In order to study the stability of such waves, we actually consider the much more general 
framework of stability/instability of operator pencils of the form 

L(A) = X 2 Id + 2c\d x + H, 

where H is a fixed self-adjoint operator, acting on L 2 , which satisfies certain mild constraints. 
This question has been considered before, mainly in a pure operator-theoretic sense. Note the 
work Chugunova-Pelinovsky, [8] in which they consider similar problems, but where the operator 
H is a finite symmetric matrix. In our work, we will need the general theory developed in the 
works of Shkalikov, [33], see also Azizov-Iokhvidov, [3] for a thorough review of the available 
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results. We also develop some of the theory below, as we shall need it for our purposes, see 
Section 3.1 below. 

Going back to the problem at hand, we consider linear, second-order in time equations in the 
general form 

(18) utt + 2uu tx + Hu = 0, (t, x) £ R 1 x R 1 or R 1 x [-L, L] 

where H = H c is a self-adjoint operator acting on L 2 , with domain D(H) and w is a real 
parameter. Note that it is better at this point to consider uj as an independent parameter, and 
to ignore the fact that in the applications to = c. 

Based on our examples, we saw that it is reasonable to assume that the self-adjoint operator 
H has one simple negative eigenvalue, a simple eigenvalue at zero (which is naturally generated 
by (p') and the rest of the spectrum is contained in (0, oo), with a spectral gap. Even though our 
approach is quite systematic and should be able to handle a more general situation 6 , we choose 
to implement these assumptions for simplicity of the exposition. Thus, we require 

' a{H) = {-<5 2 }U{0} Ua + (H), a + (H) C (a 2 ,oo),a > 
r , I H<f> = -5 2 <f), dim[Ker(H + 5 2 )] = 1 

^ ' | P^o = 0, dim[Per(P)] = 1 

, U\\ = Ho\\ = i 

We shall need the following spectral projection operators 

P :L 2 ->{0}- L ;P o /i :=h-{h,4>)(t> 

P 1 :L 2 -> {</>, ^o^; Pih :=h- (h, 4>) - (h, V> } V'o 

Our next assumption is essentially that Hi is of order higher than one. We put it in the following 
form: for all r >> 1, we require 

(20) (H + r)-^ 2 d x (H + r)- 1 / 2 , (H + r)" 1 ^ G B(L 2 ), 

Note that the quantities in (20) are well-defined, since for all r >> 1, we have that H + r > 
(r - 5 2 )Id > 0. 

An easy consequence of (20) is that H~ 1 P\d x P\ G B{L 2 ) as well. This follows easily from the 
resolvent identity, since H^P^Px = Pi(H + T)~ l d x Pi + tH^P^H + T) -1 d x Pi. In addition, 
the following non-degeneracy condition is also required 

(21) (0',Vo)^O. 

Lastly, we assume that H has real coefficients. We formulate as follows 

(22) ~Hh = Hh. 

An important observation, that we would like to make right away (and which will be used 
repeatedly in our arguments later on) is that for every A > 0, the operator (P+A 2 ) : {(f)}~ L — > {(/)}~ L 
is invertible, since (H + A 2 )||^}± > \ 2 Id. 

The following theorem is the main result of this paper. 



6 for example multiple eigenvalues at zero, which is the case in higher dimension 
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Theorem 1. Let H be a self-adjoint operator on I? . Assume that it satisfies the structural 
assumption (19), (20) as well as the non- degeneracy assumption (21) and (22). 

Then, if <^i? — 1 [V^o] ? V'o) — ®> we have instability (in the sense of Definition 1) for all values of 

Otherwise, supposing (ff _1 [^o] 5 V'o) < 0; we have 

• the problem (18) is unstable if uj satisfies the inequality 

(23) 0<U< — 1 ==:u*(H) 

• i/ie problem (18) is stable, if uj satisfies the reverse inequality 

(24) |w| > 
Remarks: 

• Note that the result is a complete characterization of the stability and instability prop- 
erties of the abstract quadratic pencil problem (18). In essence, it says that there is a 
critical number ui*(H), below which the problem is unstable and above which, there is 
stability. 

• Note that the critical number ui*{H) may also depend on 7 c, so that the equality c = 
u*(H,c) (as in (17)) will be an implicit one. In any case, the conclusions (23) and (24) 
hold true, and should be used to determine stability, regardless of this implicit dependence. 

• A more succinct way of defining ui*(H) is the following 

r +oo if <#$>,$i> >o 

u*(H) = l ^ if <#^,y^)<o 

{ 2^/-<H-i[^W,',> 

Then, we can characterize stability in terms of the inequality as follows: \uj\ > oj*(H). 

An immediate corollary is the following "spectral stability/instability" type of statement. We 
can write (18) in the form 

u ) =( 1 )( u )=-r( u 

u t J t \ -H -2ud x J \u t J \ u t 

Corollary 1. In the statement of Theorem 1, assume in addition that [Hh( — )](x) = (Hh)(—x). 
Then, in the cases of instability, there is X > 0, so that A, — A are both eigenvalues of T and 
moreover 

a(T) C {AjUj-AjUiR 1 . 
// on the other hand, there is stability, we have a(T) C iH 1 . 

2.2. Applications. 

2.2.1. Stability for traveling waves of the "good" Boussinesq equation. The linear stability/instability,! 

for the case p = 2, was considered in the literature. The linear instability for |c| < was 
shown by Fal'kovich-Spector-Turitsyn, [12]. A more complete analysis, which handled both the 
stability and the instability regimes was performed by Alexander and Sachs, [2]. The idea behind 
the proof is based on an Evans function calculations, but the authors needed some computer 
assistance 8 . 



7 and in fact, for the applications in mind, H will have the form of (14), where the dependence on c is pretty 
transparent 

8 We were also informed, [1], that Alexander, Sachs and Pego had a different proof, which is unpublished. 
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Theorem 2. The traveling wave ip c of the Boussinesq equation (10) is linearly unstable, ifp > 5. 
-tf 2 < p < 5, then it is linearly unstable ifO < \c\ < ^ p 2 1 and linearly stable, when ^ p 2 1 < |c| < 1. 

Remarks: Note that the result in Theorem 2 matches precisely the results for orbital stability 
for the "good" Boussinesq model. Recall that orbital stability was shown by Bona-Sachs for 
|c| > ^ P 2 1 an d the nonlinear instability by Liu, [25]. 

The critical case |c| = ^ P 2 1 provides a (marginally) linearly stable situation, according to our 
result, but note that Liu, [25] shows nonlinear instability there. This is due to the fact that in 
this particular case, there is an additional eigenvalue at zero (which is unaccounted for in terms 
of symmetries), which is responsible for a secular nonlinear instability. 

2.2.2. Stability for the Klein- Gordon- Zakharov system. 

Theorem 3. Let c G (—1, 1). Then, the traveling wave solution (ip(x — ct),ip(x — ct)) described in 
(9) is spectrally /linearly stable for \c\ £ p^, 1) and linearly /spectrally unstable for \c\ G [0, ^). 

Remarks: Note that the linear stability results match precisely the orbital stability results 
in [7], except at the endpoints |c| = At this point, we have linear stability, according to 
Theorem 3, but it is unclear whether the wave is orbitally stable or not. 

2.2.3. Stability for the beam equation. We state the relevant results for the beam equation. In 
it, we need to make some assumptions regarding the existence of such waves and the spectrum 
of the corresponding operator H, defined in (14). We should note that the existence of ground 
states, in a sense to be made precise below, was proven by Levandosky, [26]. Note however 
that our results apply to both the periodic and the whole line cases, with appropriate boundary 
conditions in each case. 

Theorem 4. Let p > 3 be an odd integer and L C (— \/2) be an open interval. Assume that 
there exists a one parameter family {(p c }cei of solutions to (12) (where L>0 could be finite so 
that 

• ip c £ H 1 and c — > \\^>' c \\l 2 is a differentiable function on I. 

• The operator Tl c = d% + c 2 c^ + 1 — ptp p satisfies (19) and (21). 
Then, the wave (p c is linearly stable if and only if d c \\(p' c \\ < and 



Ml/-" 

c > 



\<Pc\ 



-2WJ" 

That is, if d c \\(p' c \\ > or d c \ |v?' c | | < 0, but \c\ < _ 2 g P fj|j,/ \\ > the wave (p c is unstable. 
Remarks: 

(1) Levandosky, [26] has shown (in the case of the whole line) that a family of "ground state" 
solutions may be constructed via a variational procedure as follows. More precisely, one 
first solves the variational problem 

J ( Z ) = I-oo( Z lx - c2 4 + z2 ) dx -> min 

subject to I(z) = z p+1 (x)dx = 1. 



(25) 



i 



Then, if one sets <p c := J(z c )p- 1 z c , where z c is the solution to the minimization problem, 
then ip c solves precisely (12). Note here that the restriction |c| < y/2 assures that such a 
problem has a solution and in fact J(z) > 0. Similar approach works in the periodic case 
as well, [11]. 
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(2) It should be noted however, that the ground state solutions produced via (25) need not be 
the only way of producing the family {(f c } and in fact it is not, [11]. One reason for that is 
the apparent lack of uniqueness in solving this type of minimization problem. Second, one 
may have "excited" states (or local minima of (25)), which are also interesting solutions 
to consider. We have not checked out their stability, but it is likely that some of them are 
stable. We also point out that the variational methods for studying the (orbital) stability 
of ground states are inapplicable for these excited states, whereas Theorem 4 gives an 
exhaustive answer. 

(3) We have shown in Section 1.2.1 that all reasonable families of solutions must have the 
property that % c has a negative eigenvalue and a zero eigenvalue, with the eigenvector 
tp' c . The simplicity of such eigenvalues is however difficult to prove theoretically, due in 
particular to the lack of explicit formulas for ip c . Numerically however, one can easily 
check that T~L C has simple negative and simple zero eigenvalue, even for most of excited 
states, [11]. In fact, in all numerical runs in [11], 7~L C never failed to satisfy (19) and (21). 

(4) The restriction that p is odd is technical and it is likely to be unnecessary. More precisely, 
the issue is only the differentiability of the function x — > | </?(#) | p_1 <£>(x). 

The paper is organized as follows. In Section 3.1, we first present the relevant basic results of the 
Shkalikov's theory for quadratic pencils, after which we introduce our main tool, the function Q, 
list its main properties and derive the proof of our main result, Theorem 1. This is done modulo 
the somewhat technical Lemma 2 (which establishes the linear stability at w = u*(H)), whose 
proof is given in Section ??. To achieve this, we study and establish the Laurent expansion of 
the function Q at zero, which may be of independent use in future investigations. 

In the remaining sections, we apply the main result to obtain sharp results for linear stabil- 
ity/instability to the various examples announced above. In Section 4 we consider the Boussinesq 
model and prove Theorem 2. In Section 5, we set up the Klein-Gordon-Zakharov system and find 
the set of speeds that yields stable traveling waves as stated in Theorem 3. Finally, in Section 6, 
we establish Theorem 4 regarding the stability of the traveling waves for the beam equation. 

Acknowledgements: We would like to take the opportunity to thank our colleagues and 
collaborators Sevdzhan Hakkaev, Panos Kevrekidis and Dmitry Pelinovsky for the numerous 
discussion on related topics. 

3. Proof of Theorem 1 

We first start with a preliminary material, which is the Shkalikov theory for quadratic pencils, 
[33]. 

3.1. Some aspects of the Shkalikov's theory for quadratic pencils. We follow the pre- 
sentation of Shkalikov, [33]. We are especially interested in his index formula, which relates the 
number of unstable eigenvalues of H with the number of unstable eigenvalues of the pencil defined 
in (17). 

Following Shkalikov, [33], introduce a quadratic operator pencil in the form 

.4(A) = \ 2 F + (D + iG)X + T, 

where the coefficients F,D,G,T are operators on a Hilbert space H, satisfying the following 
conditions 

(i) F - bounded, invertible and self-adjoint 

(ii) (T, Dom(T)) - self-adjoint, invertible 

(iii) D > 0,G are symmetric; Dom(D), Dom(G) C Dom(T) and D,G are T— bounded oper- 
ators. 
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We say that an operator M is T— bounded, if 

irr^Mirr 1 / 2 G 13(H) 

We say that £ £ /o(^4), if ^(0> w ^ n domain Dom(T) is invertible. 

Clearly, in our application, we will be interested in the case F = Id, D = 0, G = —iojd x , 
T = H. However, note that T = H is not invertible in our case of interest. Nevertheless, 
we introduce the associated quadratic pencil 

i(A) := X 2 F + X(D + iG) + J 

where F = ^-^F^' 1 / 2 , G = \T\- l / 2 G\T\- 1 / 2 . We will be also interested in the spectrum 
with respect to a smoother space. Namely, introduce the space with a norm 

\\x\U_, := IIITI^Hh 

It is shown in [33] that the spectrum of A in H-i coincides with the spectrum of A considered 
on the space H. Note that in Definition 1, we only consider smooth enough solutions ip anyway. 
Thus, we need to count the number of "unstable" eigenvalues of A. 

The main result in the work of Shkalikov is Theorem 3.7. The following statement is a corollary 
of it. Here we have just presented a weaker version of the result, which will suffice for our purposes. 

Theorem 5. (Theorem 3.7, [33]) Suppose the coefficients of the pencil A, F,D,G,T satisfy 
conditions (i) — (iii) listed above. Let the numbers of negative eigenvalues of F andT, v(F),v(T) 
respectively is finite. 

Then, the spectrum of A(\) in the open right-half plane C r = {z : $lz > 0}, considered upon 
the space H^\ consists of normal eigenvalues only. Moreover, the total algebraic multiplicity of 
all eigenvalues lying in C r satisfies 

k(A) < v(T) + v(F). 

Recall that in our case however, the operator T = H is not invertible. This case is also covered 
by Shkalikov, see Theorem 4.2, [33] under the structural assumption (19). Our proof below just 
retraces his argument. 

Indeed, since Ker(H) / {0}, one needs to consider H T := H + rid for < r << 1, so that 
Ker(H T ) = {0}. 

In order to apply Theorem 5, we need to check that F = Id,G = —2iojd x are H T bounded. 
This amounts to showing that \H T \~ l l 2 d x \H T \~ l l 2 G B(L 2 ). This is a direct consequence of (20). 
We can now apply Theorem 5 to A T to conclude 

k(A T ) < u(H T ) + v(Id) = 1, 

for all small enough r > 0. Since the eigenvalues depend continuously on r (see [21], Chapter 7), 
we take a limit as r — > 0+ to get the desired inequality k(A) < 1. 

To recapitulate, we have shown that for fixed real uj and under (19), the equation 

(26) AV + 2Aw^' + Hip = 0, (A, tp) G H\ x L 2 , 

has at most one solution with KA > 0. 

Finally, it is easy to see that in this case, the solution A, if it exists, must be real. Indeed, 
suppose A 7^ A is an eigenvalue for the pencil. That is, there is tp 6 D(H) : X 2 tp + 2Xujip' +Htp = 0. 
Taking a complex conjugate and taking into account (22), we see that 

\ 2 vl) + 2\ujvJ}' + Htp = 

Thus, (A,V>) is another solution to (26),with KA > and hence k(A) > 2, in contradiction with 
the inequality k(A) < 1. Thus, A must be real and we have established 
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Corollary 2. For u G R 1 , the equation X 2 ip + 2ujXip' + Hip = 0, X £ C,ip G D{H) has at most 
one solution (X,ip) with !RA > 0. Moreover, such a pair will have X real, A > 0. 

3.2. Detecting instabilities: Proof of (23). Clearly, if uj = 0, then (18) is unstable, in the 
sense of Definition 1, as in this case A = S > and ip = <p will be a solution to (17). So, assume 
henceforth that 

We start with an useful Proposition, which appears as Theorem in [5], but it may in fact be a 
well-known result. 

Proposition 1. Assume that A is a closed, densely defined (not necessarily self-adjoint) operator 
on a Hilbert space, which is bounded from below (in.f u eD(A):\\u\\=i {Au,u) > — oo). Define its self- 
adjoint part H = 3L4 = \[A + A*] . Then 

inf{KA : A G a(A)} > mf a(H). 

In particular, if H > 0, then A is invertible. 

In (18), we use the ansatz u(t, x) = e xt ip(x), whence one gets the following equation for ip 

(27) X 2 ip + 2wAV>' + Hip = 

Our first observation is that any nontrivial solution to (27) will have a non-trivial projection onto 
cp. Indeed, assuming that ip _L <p and taking into account that {(p} 1 - is invariant under H, we can 
apply the projection operator Pq on (27). We get 

(28) (H + A 2 + 2LoXP d x P )ip = 

Considering now (H + A 2 + 2ujXP d x P ) : {<p} L -)• {<p} L , we see that (H + A 2 )| w ± > A 2 , 
whereas 2uXPod x Po is skew-adjoint operator on {<p}' L . It follows that %l((H + X 2 + 2uXPod x Po) = 
(H + A 2 )||^}± > X 2 ID, whence (H + A 2 + 2uXPod x Po) is invertible, by Proposition 1. It follows 
from (28) that ip = 0. Thus, (ip,</>) ^ for any non-trivial solution of (27). 

Since solutions of (27) are up to multiplicative constant, we may take ip : (ip,<p) = 1. That is 
ip := (p + v , where v _L <f>. We get the following equation for v 

(29) (A 2 + 2uXd x + H)v = (5 2 - X 2 )(p - 2uX(p' 

Taking dot product with (p (and taking into account (<p' , (p) = and ((H + X 2 )v , (p) = 0) yields 

2ujX(v',<p) = 5 2 - X 2 

\2 r2 

<"■*> -^f 

Taking Po i n (29) (and observing that Po u = u ) ° n th e other hand implies 

[H + A 2 + 2wAP §rPo]w = -2uX<f/. 

We now observe that the operator C\ = H + X 2 + 2ujXPod x Po : {(p} 1 - — > {(p} 1 ' is invertible by 
Proposition 1 and the representation 

C X = (H + A 2 ) + skewsymmetric : {cp} 1 -> P + A 2 | {0} ± > A 2 /d > 0. 

Thus, v = -2ujX[H + A 2 + 2uXP d x P }- 1 [(p'} G {^}- L , since 0' G {0}^. 

From this analysis, we can say that the spectral problem (27) has a real solution A if and only if 

(30) 0(w;A) := ([H + X 2 + 2uXP d x P }- 1 [<P'], ( p') + ^r^- 
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has a positive root 9 A > 0. In view of Corollary 2, we have that this condition is necessary and 
sufficient for instability of (18). We have proved the following 

Proposition 2. If H satisfies the condition of theorem 1 and Q is as in (30), then a necessary 
and sufficient condition for instability is that the function Q vanishes for some Ao > 0. 

Our next result claims the continuity of the function Q(co, A) 

Lemma 1. Let {(u n ,X n )} 6 R+ x R+ ; so that (uj n ,X n ) — > (ujq,Xq) 6 R+ x R+. Then, 

lim£?(o; n ; A n ) = ^(w ;A ) 

n 

We provide the straightforward, but somewhat technical proof of Lemma 1 in the Appendix. 
We now continue with a verification of the fact that Q(cu,-) changes sign in (0,oo) (and hence 
vanishes somewhere there), provided either — 5 2 (//~ 1 [^q], ^q) = | (<j),ip' Q ) \ 2 — 5 2 \\H~ 1 / 2 Pq[^' ]\\ 2 < 
or < \u\ < oj*(H). 

3.2.1. Analysis close to A = oo. We will show that under appropriate assumptions on H, we have 

(31) lim g{\) = — > 0. 

Indeed, lim A ^oo = ^jz, so it remains lim^oo ([H + A 2 + 2u]\P d x P ]~ 1 [</>'], (f>') = 0. 

This requires a bit of elementary analysis of the operator C defined above, in the limit of 
A — > oo. We have already established that C\ is invertible on {(f)}- 1 . We now need an estimate 
on the norm of its inverse. 

Proposition 3. For every A > 0, we have H + A 2 + IujXPqBxPq : {<f\' L — > This operator 

has an inverse (in the said co-dimension one subspace) and its inverse satisfies the estimate 

(32) \\(H + A 2 + 2u\P d x P )- 1 \\ w ^ m ± < \- 2 

Proof. We have already checked the invertibility. Let g G {4'}' L is real-valued arbitrary function 
and / = (H + A 2 + 2oj XPodxPo) -1 ^] G {(p} ± , (note / is real valued as well), so that 

(H + A 2 + 2u\P d x P )f = g, 

Taking dot product with / yields 

A 2 ||./|| 2 < ((H + X 2 )f,f) = ((H + X 2 + 2u;XP d x P )f,f) = (f,g) < \\f\\\\g\\, 
whence \\(H + A 2 + 2ujXP d x P )~ 1 g\\ = \\f\\ < X- 2 \\g\\, as stated. □ 

Going back to the formula for C^ 1 (which now makes sense), we observe that 
\\uvsuv\{[H + X 2 + 2uXP d x P ]- l [(t)']A')\ < ||0'|| 2 limsupA- 2 = 

whence we have established (31). 



9 Often times, we will omit the dependence of Q on to, since cj will be mostly fixed 



14 MILENA STANISLAVOVA AND ATANAS STEFANOV 

3.2.2. Analysis close to A = 0. Regarding the behavior close to A = e ~ 0, we observe first that 

e 2 - 5 2 5 2 



+ 0(1), 



4co 2 £ 2 Au 2 s 2 
We will show that 

(33) {[h + a 2 + 2„a w„ r vi. r> = h i^M^W + 0(E_,) ' 

Before we prove (33), let us pause for the moment to show that it implies all our statements for 
instabilities. First, if 5~ 2 \ (4>,ip' ) | 2 - WH' 1 / 2 P [ip' ] \\ 2 = - (iP" 1 ^], Vo) < 0, then the coefficient 
in front of e~ 2 for Q(e) becomes 

WMf S 2 4^((0, - g^V^hOH - 5 2 

l + 4w 2 ||iT-V2p [^]||2 4w 2 4w 2 (l + 4w 2 ||F- 1 /2p [^]||2) 

whence lim e ^o+ G{s) = — oo and hence Q changes sign in (0, oo) and it has a root Ao > 0, hence 
instability. 

If, on the other hand one has|(0,V o }| 2 ~ ^\\H~ 1/2 PqWqW = ~* 2 (H^Wol^'o) > 0, the 
solution to the inequality 

W,1>o) 2 ^ 2 c0 

l + 4w2||F- 1 /2p [^]||2 4w 2 

in terms of oj reads as < \oj\ < oj*(H), that is (23). But we observed that this implies change 
of sign for Q. Thus there is Ao > 0, so that Q{\q) = 0, hence instability. 
Back to the proof of (33), we need to compute the leading order term in 

((H + e 2 + 2uePod x P y l [(/>'}, <j>') . 

Let z G {4>} ± solves (H + e 2 + 2ujePod x Po)z = 4>'. Decompose z = aipo + q, q € {(f), 4>o} ± - Observe 
that q = P±q. The equation for z becomes 

(34) ae?ifo + 2aweP [if/ ] + (H + e 2 + 2ueP d x P )q = <j>' 

The equation again has two different components - along ?/>o and perpendicular to it. 
Along the direction of tpo, we take dot product with ipQ to obtain 

as 2 + ({H + e 2 + 2ueP d x Po)q,ipo) = {tf,ifo) 
which implies (note (q,ipo) = 0, since q = P±q _L i/jq) 

(35) as 2 -2ooe(q,^ ) = (<jf ,-fo) 

Along the orthogonal direction, we take Pi in (34) - note that P\Pq = P\ and P\Po[ip' ] = 
P Q P 1 [iP' Q } = PoW ]). We get 

(H + e 2 + 2ueP x d x P{)q = Pi[(jJ] - 2aujeP [^] 

Denoting the self-adjoint operator Hi = H\Pi = P\HP\ > a 2 and consider C £ = Hi + e 2 + 
2ujePid x Pi. Since Hi + e 2 > a 2 , we may write 

C £ = Hi[Id + H~\e 2 + 2uePid x Pi)], 

whence by (20), lim e ^o+ \\Hi 1 (e 2 + ^ u,£ Pi9 x Pi)\\ l 2 ^l 2 = 0) we can construct C^ 1 in terms of 
Neumann series. In fact 

C^ = H^ + 0(e). 

Since H~ 1 = O(l), it follows that 

(36) q = £- 1 [-2aweP [V'o] + p iW\\ = ~2aueH^ P [$,] + 0(1) + aeO(e). 
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Expressing q from (36) back in (35) yields 

ae 2 - 2uje{-2auje (H^P Q [^ Q ], ij/ Q )) + 0(e)(l + ae 2 ) = (</>', Vo) 

Note 

(flr^M,^) = (h^pm,p w }) = \\H~ 1/2 [PM]\\h- 

As a consequence, 

ae 2 (l+4a; 2 || J ff 1 - 1/2 [Po[V'o]]|l! 2 +0(^) = (^,V'o> + 0(e), 

whence 

4 .<?•*■> + o(«-) 



e 2 i _i_ /i,.,2ii it -1 / 2 rp_r.// 11112 



l + 4a; 2 || J ffr 1 ^[Po[^]]|li 2 
Going back to (36), we conclude that q = 0(£~ r ). Now 

((H + e 2 + 2ujeP d x Por 1 4>'A') = a (Vo, <//) + (<?,</>') = 

1 <</>', V>o} 2 



^1 + 4^11^— /- [jPo[ ^]]||2 2 

which establishes (33). 

We now progress to the proof of (24). Introduce the notation 

junstable ._ ^ . (lg) j g ungtable } 
A stable ._ | w . ^ ig gtable } 



+ 0(e- 1 ). 



We have already shown that A unstable / 0, in fact (-u*(H),u*(H)) C A unstable . 

The rest of the argument proceeds in two steps. First, we show that (u*(H), oo) C A stable and 
then, we prove that u*(H) € A stahle . 

3.3. Proof of (-oo,u*{H)) U (w*(fT),oo) C „4 staWe . Let w : > w*(fT) > 0. As we saw, this 
condition implies that the function Q{ujq; A), defined in (30), has 

lim Q(ujq;\) = > 0, lim G(wq;\) = +oo. 

A^oo 4w 2 A^0+ 

By the continuity of the function A — > Q(ujq;\), there are three options. The first one is that 
Q{ujq]\) > for all A > 0. By Proposition 2, this would mean that Q(uq,X) never vanishes and 
hence, (27) does not have solutions, whence stability for ujq follows. We claim that this is the 
only viable option, by disproving the other two. 

The second option is that Q{ojq\ A) has at least two different roots in (0, oo), say < Ai < A2 < 
00. This implies that (27) will have two different solutions, (Ai,V>i), ^2,^2), thus contradicting 
Corollary 2 of the Shkalikov index theory. We can therefore refute this option. 

The third possibility is that G{ojq; A) has a double root, say Ao = A (^o) > 0) brrt otherwise is 
a non-negative function. If we manage to rule out this possibility, we will be done. This in fact 
not very hard. Assume for a contradiction that the double root actually occurs at Ao- Let ji 7^ 
and consider the associated perturbed problem for := H — [i 2 (•, <j)) 4>. Clearly, PqH^Pq = H 
and thus, according to (30) 

\2 x2 2 2 

(37) £tf>, A) = ([PqH^Pq + A 2 + 2LoXP d x P ]- 1 W], 4>') + ~^~ 2 - = S(u, A) - 
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The analysis that we performed for H remains valid for H^, at least for < fi « 1. That is, the 
root (if any) of A — > Qh A) will give exactly the instability eigenvalues of H^. In particular, 
by the Shkalikov's theory, the number of zeros of GhAw, A) is either one or zero. 

But if Q{ujq, A) has double zero at Ao, it follows by (37) that for all \i small enough, the function 
A — > Vh^^loo, A) will have at least two separate zeros Ai < Ao < A2, A2 — Ai = 0(p?). In any case, 
as we have alluded to above, that is a contradiction with the Shkalikov's theory, which requires 
there is at most one unstable eigenvalue for H^. With that, we have reached a contradiction, 
which implies that double roots are impossible for A — >■ Q(oj, A), when uj > 00* (H). Hence, 

(w*(fl"),oo) C A stable . 

3.4. Proof of u*(H) £ A stable . 

Lemma 2. The problem (18) is stable at uj*(H), i.e. oj*(H) G A stable . 

We now proceed to show Lemma 2. Denote in this section loq := oj*(H) for sake of brevity. 
We shall need several propositions. 

Proposition 4. For any uj, the function e — > e 2 G(ui,e) is real analytic (close to zero) and it has 
the Laurent expansion 



(38) g(u,e) =e- 2 D_ 2 {uj) + ^2D 

where D-\{uj) = 0, 



3 



£>_ 2 (w) 



3=0 



l + 4w 2 ||i7- 1 /2p [^]||2 4u; 2' 

and the functions {Dj}'*L l are smooth functions of uj. Finally, the radius of analyticity r{ui) 
satisfies 

inf r(uj) > r ab > 0. 

oje[a,b]c(0,oo) 

That is, on every compact interval I = [a.b], one may choose a common radius of analyticity. 

Proof. Most of the statements in Proposition 4 were in fact considered before, see for example 
the formula for D_2(u;), derived in (33). 

Let us first indicate the real analyticity of the function e 2 Q(e). We have, from the defining 
equations and in our previous notations, 

1 5 2 

G{uj, e) = a £ <Vo, <!>') + (q e , 4>') + ^ ~ 4^2 ■ 
where a = a £ and q = q £ are determined from the pair of equations 

(39) as 2 -2uje(q,^' ) = (<f/, ^ > 

(40) q £ = C-\-2aeu>Pv[i>' G ] + P\[4>'\), 

where C £ = (H + e 2 + 2ujePid x P\). As we observed before, by the von Neumann expansion, 

00 

C- 1 = J2(-z) k (H-\e + 2ujP 1 d x P 1 )) k H^ 1 = tff 1 + 0(e) 

k=0 

and is real-analytic in e in a region |e| < £0 << 1- Substituting the expression for q in the 
expression for a yields the formula 

as 2 + We 2 {C- l P W Q ], [$,]) = 2UJ8 (C^P^], [^]> + Vo> 
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and 



ae 



2 



l + 4^(£-ip [^],[^]> 

Observe that <(£~ 1 Po[V'o]' V'o) = ll-^ -1 ^ 2 -^^] II 2 + 0(e), whence the denominator stays away 
from zero, whence the real analyticity of e — > e 2 a(e). Furthermore, the coefficients Dj(ta) (which 
can be written by means of the Cauchy formula) are smooth functions of uj. Substituting now 
this back in the formula (40) yields that e — > eq(e) is real analytic too. Observe also that the 
radius of analyticity of both functions may be estimated below in terms of Therefore, it is 
uniformly bounded away from zero, when to belongs to some compact interval, away from zero. 
We now compute D^i^), D_i(uj). To that end, introduce A, B, so that 

a(e) =£- 2 A + £- l B + 0{l). 

We compute q up to order 0(1). We have C~ l = H^ 1 + 0(e) and therefore 

q = (H' 1 + O(e))(-2e~ l Au;P [^ ] - 2BuP [i> ] + P^']) = 

= -2^we- 1 F- 1 P [V'o] - 2BujH- 1 P [i)' Q ] + H^P^'} + 0(1) 

Plugging this in (39), we get two equations - at the scale of e° and e 1 

A - 2u(-2Au (H^Poty'o},^)) = (<f>', Y>o) 

B - 2u((H- 1 P l [<t ) , )^ , Q ) - 2Bu (H-iPoM],^)) = 0. 

From those, we get 

<</>', rk) 



A 



l + 4a; 2 ||P- 1 /2p [^]||2 

2a;(g- 1 Pi[0 / ],^ o ) 
l + 4 W 2 ||i7- 1 /2p [^]||2 

We are now ready to compute Q(oj,e) up to order 0(1). We have 

1 S 2 

G(oj,e) = a e {4>\ilJo) + {q £ A') + -^ 



4lo 2 4u> 2 e 2 

2 



= e- 2 (a ((f/, iPo) - ) + e- 1 (B (<//, ^) - 2Aw {H~ l P [^ 4>')) + 0(1) 

It remains to observe that 

D_i(w) = B ((/)', 4i )-2Au(H- 1 P [<iP' }, ( t } ') =0. 

□ 

Several things to note. First, recall that the definition of ujq = co*(H) is so that Z?_2(wo) = 0. 
In fact, D_2(<^) < 0,lu < ujq and D_2(^) > 0, w > ojq. 
Now, consider the real-analytic function 

00 00 
G(uo,e) = e- 2 D_ 2 (ou ) + J2 D J^o)e j = J2 D j("o)e J . 

3=0 j=0 

By the analyticity, we have that there is k, so that D^loo) 7^ 0. Indeed, otherwise, the real- 
analytic function e — > Q(coo,e) will be identical to zero, at least for all e small enough. This is 
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however a contradiction with the meaning of the function Q(ujq,-). Indeed, not only it implies 
instability at ujq, but it gives eigenvalues A = e G (0, So), clearly in contradiction with our earlier 
conclusion (basically a consequence of Shkalikov theory) that there is at most one eigenvalue. 
Thus, there is k : D^q) / 0. Let ko = min{/c > : D^ujq) / 0}. 
Claim: Dk (ujo) > 0. 

Indeed, take ojj — > ojq+. Since ojj G A stable , take any e : < e < ro = inf r(uij) (which is strictly 
positive by Proposition 4), 

fco-i 

< limsu P g(o; J -,£) = limsup^- 2 !^^) + (^ D m {uj)e?) + D ko ( Uj )e k ° + O(e ko+1 )) = 

3 3 m=0 

= D ko (u )e ko + O(e ko+1 ) 

Clearly, the last inequality will be contradictory for small e > 0, unless Dk Q (u}o) > 0, which was 
the Claim (recall Dk {ooo) 7^ 0). 

Having established this claim, we easily finish the proof of our claim that u G A stable . Indeed, 
we have lim^oo ^(^o, A) = j^j > and lim^o+ ^~ k °G(uo, A) = D ko (u)o) > 0. Again, assuming 

the contrary (i.e. uoq G j[ unstable ^ implies the existence of Ao > 0, so that Q{uiq, Ao) = 0. By the 
behavior of A — > Q{ojq, A) at A = 0, oo and the fact that the function G(ujq, A) cannot have more 
than one zero in (0, oo), we again conclude that QiuiQ,-) has double zero at Ao- As we saw in 
Section 3.3, this leads to a contradiction, if we consider the function A — > Gh^Iloq, \). Namely, 
for all sufficiently small /x, the function Qh^q-.X) will have at least two zeros, in contradiction 
with the Shkalikov's theory. 

3.5. Proof of Corollary 1. In the course of our exposition for the Shkalikov's theory, we have 
proved that the spectrum of the pencil (and hence the spectrum of the operator T) is invariant 
under complex conjugation. More specifically, if solves (18), then so does (A, ip). 

Now, under the assumptions in Corollary 1, we can also show that the spectrum is invariant 
under the transformation A — > —A. Indeed, if (A, ip) solves (18), then so does (— ■))■ 

Thus, in the unstable situation, we have shown that there is unique positive eigenvalue A of T 
and hence there is an unique negative one as well. Moreover, this argument excludes the possibility 
for more negative spectrum, since this would imply more positive (unstable) eigenvalues for the 
pencil, which is forbidden by the Shkalikov's theory. 

In the stable situation, one shows by contradiction argument, similar to the one in the previous 
paragraph, that there is no spectrum off the imaginary axes and thus <x(T) C zR 1 . 

4. Linear stability analysis for the Boussinesq model 

The proof of Theorem 2 follows directly from Theorem 1, applied to the linearized problem 
zu + 2cztx + Hz = 0, where H is as defined in (16). In order to follow this program, we 
need to verify that the conditions in Theorem 1 for the operator H are met, then we need to 
compute the critical index ui* (H) and finally, we need to establish that the linear instability for 
zu + tcztx + Hz = is equivalent to the linear instability for the actual problem (15). We will 
do this in a sequence of several propositions. 

Proposition 5. (Spectral properties of H) The operator H, defined in (16) satisfies (19), (20), 
(21), (22). 

Next, 



Proposition 6. The critical index oj*{H) = ^^gSi ^ 
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Regarding the existence of the unstable solutions, we have already explained that an instability 
in Zu + tcztx + Hz = implies instability in (15). We need to show then the reverse. 

Proposition 7. Suppose that V £ il^R 1 ) solves the equation 

X 2 V + 2c\V + T[V] = 0, 
for some A > 0. Then there exists a function Z G H^^R 1 ), so that Z' = V and so that Z solves 

X 2 Z + 2c\Z' + H[Z] = 

Let us show now that Theorem 2 is a simple consequence of Propositions 5, 6 and 7. Indeed, 
since we can apply Theorem 1, we only need to find the intervals, in which the speeds yield stable 
traveling waves. To that end, we set up the inequality 



V 5 -p 

The solution to this inequality is 

l>|c|>^H. 

In the complementary set, < |c| < we have instability. This finishes the proof of Theorem 

2, modulo the claims of the Propositions. 

4.1. Proof of Proposition 5. We start off by noting that (22) is obvious by inspection. Sim- 
ilarly, (20) follows from the fact that H has one negative eigenvalue 10 and thus, for r >> 1, we 
have that {H + r) > Id, hence is invertible and moreover 

{H + t)- 1 : L 2 — >■ H 4 , 

and more generally for all s G [0,1], {H + r)~ s : L 2 — > H As . Therefore, we conclude (H + 
T)~ l / 2 d x (H + t)^ 1 / 2 , d x (H + t) -1 : L 2 — > L 2 are a bounded operators. 

4.1.1. Proof of (19). To that end, note first that H can be written in the following form 

h = -d x cd xi c -.= -d 2 £ + (i - c 2 ) - ptfr 1 

The operator £ is a second order differential operator, which appears in the study of the KdV 
equation in exactly the same way - as a linearization around the traveling wave solution ip c . 
Its spectral properties are well-documented and we just remind the ones that are of concern to 
us. Namely it has a simple eigenvalue at zero (with eigenvector (p' c ), single and simple negative 
eigenvalue and it has a spectral gap, that is the rest of the spectrum lies in a set [k, oo), k > 0, 
see for example [6], Section 3. 

We will now be able to infer the same structure for the spectrum of H, mainly due to the 
representation H = —d x Cd x . First, it has an eigenvalue at zero, with eigenvector (p c {H(p c = 
— d x £{(p' c ] = 0). This eigenvalue is simple (otherwise a contradiction with the simplicity of the 
zero eigenvalue for £). 

Next, we check that H has a negative eigenvalue. Taking two spatial derivatives in (6) we get 

(4i) v'c" - (i - c 2 )^+p^\"+ P {p - i^r Ve) 2 = o 

taking dot product with Lp" c and taking into account that —C(p" = <p"" — (1 — c 2 )^" + pif p c ~ 1 ip" , 
we get 

- </y c ', ^> + p{p - 1) J ^- v c )y c 'dz = o 



10 to be proved below in this proposition 
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But JVrVc)V c 'dx = -2^ J>rVc) 2 dx. Thus 

Note that if p > 2, the expression is negative, thus yielding directly (by the Ritz-Reileigh criteria) 
that the minimal eigenvalue 

V>:|HI=1 

Even in the case p = 2, this argument implies the existence of a negative eigenvalue. Indeed, 
assuming that H > (and since we already know that (p c is a simple eigenvalue), it follows that 
-^I{vp c } ± > 0- Thus, since G {i^c}- 1 , it should be that (H<p' c ,<p' c ) > 0, a contradiction with 
(Htp' c , ip f c ) = 0, in the case p = 2. This shows that Xq(H) < 0. 

Finally, we need to establish that Xq(H) is simple. This would follow, if we manage to show 
that \i(H) > (and indeed, since is an eigenvalue, it would follow that \i(H) = 0). Denote 
the (smooth) eigenvector of £ by £. One may apply the Courant maxmin principle for the first 
eigenvalue, which states 

Ai(H) = sup mf 

z -LQU±z \\u\\ z 

Taking z = 4>' above yields 

wm . r (Hu,u) . (Hu,u) . (£u',u') . (£u',u') 

XAH) = sup ml > mf = mf — ——75 — = mf — — — — 

zj lqu±z \\U\\ Z u±<j>' \\U\\ Z u±<t>' \\u\\ z u'±(j> \\u\\ z 

where in the last identity, we have used that u _L 4> exactly when v! _L 4>. It remains to observe 
now that since C has a simple negative eigenvalue, with eigenvector cj>, we have £|{<^}i > and 
hence (Cu',u'} > (since u' € {(f>} ± ), whence 

A t (g)> inf ^# >0- 

Thus, property (19) is fully established. 

4.1.2. Proof of (21). Note first that in our notations, ^0 = Vc/ll^cll an d denote the negative 
eigenvector of by 0, \\<f>\\ = 1. We need to show (0', <p c ) / 0. We separate the proof in the cases 
p > 2 and p = 2. We have by our computations in Section 4.1.1 above that 

(iV c y c )<o 

It follows that 

> (H(p' c , ip' c ) = A <y c , </>} 2 + (H(<p' c - (tp' c , 4>) <j>), (ip' c - {ip' c , 4>) <f))) 

Note that (ip' c — (f' c ,4>) <p) € {4>} ± - Since Xo(H) is the only negative eigenvalue for H and it is 
simple, it follows that H\^y± > and in particular 

(h(^' c - (<p f c , 0> 0), y c - (<p f c , 0> 4>)) > 0. 

The last two inequalities imply that 

0> A <^ c ,0> 2 , 

whence (tp' c , <j>) / 0. But then (ft, ip c ) = — (ip' c , ft / 0. 
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In the case p = 2, as we have shown (H(p f c ,tp f c ) = 0, so we need to be more precise in the 
arguments above. In particular, we need to observe that since (p' c _L ip c and 4> _L ip c (as eigenvectors 
corresponding to different eigenvalues), it follows that actually 

(42) ip' c - (ip' c , 4>)4><E span[(p, ip^. 

Moreover, by (41) (in the case p = 2), we have that £[<p"] = 2((p' c ) 2 , whence 

ny c ] = -d x ny c \ = -VX' * W c - 

The last computation shows that ip' c and 4> are linearly independent and hence, we can upgrade 

(42) to 

(43) ^ ip' c - (ip' c , (p)4>e span[4>, ip^. 

By (19), in particular the spectral gap that we have established for H, and (43), it follows that 

(H(<f/ C - (<f/ c , </>} 0), y c - (<f/ c , 0> 0)> > ny c - (d, 0> 0|| 2 > o, 

where k is the size of the spectral gap. We can conclude 

Xo(<p' c ,<f>) 2 = (H^^) -(H^-^,^)^),^- (&,<{>) cf>)) = 
= - (H(<f/ C - (<p' c , <p) <P), (<p' c - (<f/ c , <p) </>)) < 0, 
whence (<p', tp c ) = - (ip' c , 4>) / 0. 

4.2. Proof of Proposition 6. Recall H = —d x Cd x . In our notations, we need to compute 

<ij-V^> = ^<#-Vc^>- 

Again, starting from (6), take a derivative in the parameter c. We obtain, 

(d c (p)" - (1 - c 2 )d c ip + pip^dcip + 2c<p = 0. 

In terms of the operator £, we have £[d c (p] = 2cip or ^d c (p = £~ 1 [(p c ]. Here recall that C is 
invertible on Ker^] 1 - = {(p' c } ± and ip c G {<p' c } ± - 

Introduce z, so that ip' c = H[z\. We know that such a z exists, since the operator H is invertible 
on the subspace KerlH] 1 - = {(p c }' L and ip' c € {<fc} ± - We have 

<ff c = H[z} = -d x C[z'} 

Thus ip c = —£[z'] + const and the constant turns out to be zero by testing this identity at x — > oo. 
Thus, z is such that £[z'] = — ip c and hence 



Now 



z' = -£ 1 [(p c ] = —^dc<P- 

1 1 

- Vo <Pc) = - ( z > Vc) = ( z 'i <Pc) = ( 9 c¥>c> Pc) = - — d c \\ip c \\ 2 . 



All in all 

/TT-l // ,/\ 1 /TT-l I /\ 1 5 C [ll^c|| 2 ] 

-{H ^) = -—(H (Pe , Ve ) = ---^- 
It remains to compute ||^ c || 2 and perform the elementary calculus operations. We have 

= m p (l — c 2 )~v~- 



2 
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Thus, 

1 9 C [||^ C || 2 ] 5-p 



4c \\ip c \\ 2 ] 4(p- 1)(1 - c 2 )' 
Clearly, for p > 5, we have (ii~ _1 V>o; V"o) = 4( p _i)"(i-c :j ) — ^ anc * nence instability. If p < 5, we get 

w vm = 1 = VCp-^i-c 2 ) 

2V-(^-Vo,^) v 7 ^ 

4.3. Proof of Proposition 7. Proposition 7 is a simple exercise. Writing down the particular 
form of the operator T, we see that V satisfy 

(d x - (1 - c 2 )d 2 x + 2c\d x + \ 2 )V = -pdl^V). 

Observe that the symbol of the operator M = d x — (1 — c 2 )<9 2 + 2c\d x + A 2 

e + (1 - c 2 )£ 2 - 2dA£ + A 2 

has a positive real part, bounded away from zero and is thus invertible on L 2 . In particular, we 
may write 

v = -phr 1 d 2 Mr 1 v}. 

Now, it suffices to set 

Z := -pM- x d x [v%- x V\, 

so that Z x = V (and hence it satisfies \ 2 Z + 2cXZ + H[Z] = 0) and check that such a function is 
well-defined and belongs 11 to L 2 (R 1 ). Since (fc~ 1 V G L 2 , it suffices to show that d x M~ x defines 
a bounded operator on L 2 or 

e 



sup 



< oo 



e 4 + (1 - c 2 )i 2 - 2d\£ + A 2 
The last inequality follows by inspection. 

5. Linear stability analysis for the Klein-Gordon-Zakharov system 

In this section, we will provide the complete details of the proof of Theorem 3. We need to first 
setup the problem at hand in the form (17). This is not straight forward, due to the requirement 
that % be self-adjoint. In fact, the choice of ^ that we have made in the formulation of (7) helps 
us accomplish exactly that. 

We take the ansatz u(t, x) = (p(x — ct) + v(t, x — ct),n(t, x) = ip(x — ct) + w(t, x — ct) for real- 
valued functions v,w. We plug it in (7) and ignore all terms 0(v 2 + w 2 ). We get the following 
linear equation for (v,w), 

v tt - 2cv tx - (1 - c 2 )v xx + v + ipv + (pw = 
w tt - 2cw tx - (1 - c 2 )w xx - (ipv) xx = 0. 



(44) 



At this point, we introduce the convenient quantity fi = \J 1 — c 2 . Setting w = z x and taking off 

-» ( v 

one derivative from the second equation, we obtain the following system for $ = 
(45) $ tt - 2S tx + n$ = 0, H:=( J ^ ) , 



11 further regularity of course holds and can be inferred in the standard bootstrap fashion 
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where 

2 

Hi = -(1 - c 2 )8 xx + 1 + V = —^ 2 d xx + 1 " ^ 

h 2 = -(i- c 2 )d xx = -n 2 d xx 

Az = ipz x ,A*v = -(<pv) x 

Clearly, the operator H is self- adjoint. It now remains to show that % satisfies the requirements 
of Theorem 1, after which, we will compute the quantity of interest (%~Vo> V'o)- We formulate 
the needed results in a series of Propositions. 

Proposition 8. The self-adjoint operator TL defined in (45) has a simple eigenvalue at zero, with 

2 

an eigenvector ('■p',—-^)- 

Proposition 9. The self-adjoint operator TL defined in (45) has one simple negative eigenvalue. 
5.1. Proof of Proposition 8. We need to solve TL ^ ^ ^ = 0. This is equivalent to 

Vf + / - + = o 

-(¥>/)' -mV = 0. 

Integrating the second equation yields g' = which we plug in the first equation. We get 

(46) - ff" + f - ^/ = 0. 



Recall now the equation (8), which defines if. In fact, taking a derivative d x in it yields 

2 m , / *V 2 , 
~H <p +<p - -^ip = 0. 

By comparing the last two formulas, we see that / = iff is a solution of (46). Moreover, by the 
standard theory for one-dimensional Hill's operators 12 (see for example [6], Section 3), iff is the 
unique solution to (46), up to a multiplicative constant. 

It remains to observe that since g' = — = we have 



2 

and thus the eigenvector, corresponding to the zero eigenvalue is (if 1 , —^)- 



1ix l 
ro < 

5.2. Proof of Proposition 9. The proof of Proposition 9 is less standard than the proof of 
Proposition 8. We set up the eigenvalue problem in the form TL ^ ^ ^ = —a 2 ^ ^ ^ for some 
a € (0, oo). That is, we need to show that there exists an unique ao > 0, so that the eigenvalue 
problem has an unique (up to a multiplicative constant) solution ( ^ ^ . Write the eigenvalue 
problem in the form 



-^/" + /-^/ + ra' = -« 2 / 

Vf - (¥>/)' = ~a 2 g. 



12 Note that this particular operator is in fact the ubiquitous L_, which appears in the linearization of standing 
waves for the cubic Schrodingcr equation, which is known to have one-dimensional kernel 
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From the second equation, we express g = d x (a 2 — ix 2 d xx )~ l [</?/]. This is possible, since 
a 2 — fi 2 d xx > a 2 Id and hence is invertible. Thus, 

9' = d xx (a 2 - At 2 <9xx) _1 [<£>/] = "4 + 4( fl2 " M^xx)- 1 ^/]- 

Plugging this in the first equation, we obtain 

V f a 2 
2 M 2J + /I 2 

To recapitulate, we have obtained the equation (47) to be equivalent to the eigenvalue problem 



(47) - fd xx f + (1 + a 2 )f - ^ / + ^(a 2 - » 2 'd^)' 1 \<p /)] = 0. 



H y ^ J = —a y g J' ^hat to P rove Proposition 9, we need to show that there exists an 
unique ao > 0, so that the operator 

M a = —fi 2 d xx + (1 + a 2 ) - ^ + ^(a 2 - 

has an eigenvalue at zero for a = ao and in addition, this eigenvalue is simple. 

Several things to note for the one-parameter family of operators M a . It is clear that M a , a > 
are self-adjoint and in addition, 

Claim: If a > b > 0, then M a > M h + (a 2 - b 2 )Id > M b . 

Let us finish the proof of Proposition 9, based on this Claim, whose proof we postpone for 
the end of this section. Denote 

A(a) = inf{A : A G a{M a )} = M (M a f, f) . 

Clearly, a — > A(a) is a continuous function and in view of Claim, a — > A(a) is an increasing 
function of its argument. 

Next, we consider the easy case a >> 1. Observe that if a 2 > sup x 3<fi 2 ^ — 1, M a becomes a 
positive operator (the operator [ip(a 2 — fi 2 d xx )~ l (f-)] is positive by its Fourier transform repre- 
sentation, see the proof of the Claim below) and thus A(a) > 0. 

The case a = presents another interesting observation. Namely, since 



M = -n 2 d xx + 1 - 



2^ 2 



we already know that Mq [(p'\ = and hence by Sturm-Liouville theory, there is a negative 
eigenvalue for Mq, that is A(0) < 0. Alternatively, one may directly compute 



1 f°° 

(M oV ?,^) = — 2 / f 3 (v)dy < o, 

M J —00 



whence again A(0) < 0. Thus, the continuous and increasing function a —> A(a) is negative at 
a = and positive for large a, whence it has exactly one zero, say ao. Thus, the eigenvalue that 
we are looking for is A = — < 0. We still need to check that this eigenvalue is simple, which 
by the equivalences that we have established means that we have to show that is an isolated 
(simple) eigenvalue of M ao . 

This is however an easy consequence of the fact that M ao > Mq + a^Id (from the Claim). 
Indeed, denote by <fo the eigenvector for Mq, which corresponds to its unique and simple negative 
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eigenvalue ([6], Section 3). Then, M |r^ i± > 0, whence by the Courant maxmin principle 
Ai(M ao ) = sup inf (M ao u,u)>al+ inf (M u,u)>al, 

z ^0 u±z:\\u\\=l u±cf> :\\u\\ = l 

Thus, we have that A(ao) = Ao(M ao ) = 0, while Ai(M ao ) > a 2 ,, which finishes the proof of 
Proposition 9. 

5.2.1. Proof of Claim. By the particular form of the operators M a , it suffices to establish that 
for all test functions /, 

(48) a 2 (Ma 2 - ^ 2 d xx )-\ V f)], /> > b 2 <[^(6 2 - ^d xx )~\^f)], f) , 

whenever a > b. This is easily seen by the PlancherePs theorem. Denoting h = iff, we have 

/°° n 2 

/°° h 2 
.oo ^ + w |ll(()|2 ^ = b2 ^ {h2 ~ ^rHMf) , 

where we have used the elementary inequality a s +4 °;^^ > ^ +4 ^^ 2 ^ , if a > 6. 

5.3. Proof of Theorem 3. Now that we have checked most of the required conditions on T-L, 

( if' \ 

let us verify the remaining ones and the quantity (T~L 1 ip' ,'tp' '). Here tpo = m I ^ , where 

m is so that HV'oll = 1- Thus, we need to compute Tl~ l I ^ ) . Recall that %~ 1 exists on 

V ~ 2 ^ 1 / 

{■^o}" 1 " an d since tpQ £ {ipo}^, h is possible to compute T-L^ 1 ^. 

This could be done almost explicitly, in any case explicit enough, so that we may compute 
( / H~ 1 ip' Q , 4>'q) precisely in terms of c. We have 

V/" + f-£f + w' = <p" 



W - (WO' = -($)'■ 



Integrating once in the second equation yields 
Plugging this in the first equation yields 



9 = — (w-2-ff)- 



(49) _, Y ' + / _V / + ^ = / 

On the other hand, recalling equation (8) in the form 



(50) -a*V + V-^ = 0, 

take a derivative in the parameter ^ in (50). We get 

q 2 3 

-aw + ^ - + = w- 
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Dividing by 2/i yields the relation 



(51) 



2/i 



+ 



2/x J 2/i 2 V 2// y /i 3 ~ ^ ' 



Comparing the equations (51) and (49) clearly implies (since we know that there is an unique 
solution by the invertibility of % on {i^o}' 1 ) that 



Thus, it follows that 



9 fi 2[ 2fi 2 



Now, 



n 



-(£)' J ' 



<p 



r' 2 V 



-2^t 2 



* 2/i 

¥>/) = 9 = 



/i 2l 2/i 2 


^ 2/i 




))■ 



/ , V / 

\ 2/i 



,<P ) + { 9 , 



<P 

2/i 2 



d,[Sy{y)\ 2 dy] d»[f<p\y)dy] , jV(j/)dj/ 



4/i 



16/i 5 



+ 



4/i 6 



Now since 



o'{y)\ 2 dy = 4 ^ {sech 1 (y / fj,)) 2 dy = 4/i /" (sech' (x)) 2 dx 
j\\y)dy = l^ jsech{xfdx, 



we have 



ft 



95 



V2^?i 



-— (^j {sech! {x)) 2 dx + ^ sech{x) 4 dx 



Recall however, that we have to also compute the normalization factor m 2 . We have 

2 
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Thus, putting the last two formulas together 



(52) 



(U- l ^^) = m 2 (u{ 



<P 



-(KY 



1 



From this identity, we get that (% _1 V>o> V'o) < f° r an values of the parameters, so in particular 
we claim that (21) is satisfied. 

Indeed, if we assume that = ((/>', ipo) = — (<j),tp' ), it would follow that ifj' € {4>} ± - Since 



H 



> 0, it follows that 



<w-y ,^o>>o, 

which is a contradiction with (52). Properties (20) and (22) are obvious and hence, we may apply 
Theorem 1. Since 



STABILITY AND INSTABILITY OF TRAVELING WAVES OF SECOND ORDER IN TIME PDE'S 27 

we have linear stability for all c satisfying 1 > |c| > /j, = \J\ — c 2 , and instability for < |c| < 
/Li = \J 1 — c 2 . Solving the inequalities yields stability for c : |c| £ [^,1) and instability for 
c : |c| G [0, ^). Theorem 3 is established. 

6. Proof of Theorem 4 

Since we cannot verify theoretically the assumptions 13 of Theorem 1, we have put the appro- 
priate assumptions, essentially requiring them, so that Theorem 1 applies. 

First, the case c = always presents instability, since H c has a negative eigenvalue. So, assume 
c > 0, the cases c < being symmetric. 

Next, we need to compute H _1 [</?"]■ This is done by just taking a derivative with respect to 
the parameter c in the defining equation (12). We have 

<?(d c <p)" + 2c^' + (d c <p)"" + (d c <p) - w?-\d c <p) = 

or equivalently H[d c ip] = —2ap". Assuming c / 0, we conclude 'H~ l Wc} = —^d c (p. Now, 

According to Theorem 1, we have instability for all c such that 

Well > 

and moreover, if d c ||</^|| < 0, we have instability for all c > 0, satisfying 

1 

c < 

or 

y/c < 



, / &=||^|| 2 



In the complementary range, we have spectral stability. 

Appendix A. Proof of Lemma 1 

The term is clearly continuous in both u,X, so we concentrate on the continuity of the 

mapping (u, A) -»• ([# + A 2 + ^APo^x-Po] -1 ^]. 

We now need to show the continuity of the map stated above. Taking a sequence (w„, A n ) — )■ 
(wo,Ao) and denoting 



we need to show that 
which follows from 



oj 

Ro = (H + \ 2 + 2oj \ Pod x P )- 1 , 



11(^-^0)^011^^2 ->0, 

which remains to be proved. By the resolvent identity, we have 

Rn - Ro = —Rn{^n ~ Aq + 2(w„A„ — W \))Po<9:cPo)Po 



13 



but, as we have alluded to before, they hold nice and steady in numerical simulations, [11] 
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and since |A 2 - A 2 , + 2(w n A„ - w A )| < C(\u n \ + \u> \ + |A | + \X n \)(\uj n - uj \ + |A n - A |), it will 
suffice to prove 

(53) limsup WRnRoPoW^) + \\R n (P d x Po)R Po\\B(L 2 ) < °°- 

n 

The first estimate follows from Proposition 3 

limsup Hi^i^ollBOL 2 ) < limsup ||i? n P ||s(L 2 ) ll-Ro-Po \\b{i?) < limsup(A~ 2 Ag 2 ) = Ag 4 . 

n n n 

For the second term, we further use the resolvent identity to write 

Rn{Pod x Po)RoPo = (H + A 2 ) 1 Pod x R()PQ — 2Lo n \ n R n Pod x Po(H + A 2 ) 1 3 x i?o-Po ) 

and we do similar expansion on the right-side with Rq = (H + A 2 ,) -1 + .... Clearly, in the above 
formula, we can pair the operators d x with various resolvents of the form (# + A 2 )" 1 , (F + A 2 ,)" 1 , 
through (20) to obtain the desired boundedness results. 
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